Fractional Brownian motion and the critical dynamics of zipping polymers 
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We consider two complementary polymer strands of length L attached by a common end monomer. 
The two strands bind through complementary monomers and at low temperatures form a double 
stranded conformation (zipping), while at high temperature they dissociate (unzipping). This is a 
simple model of DNA (or RNA) hairpin formation. Here we investigate the dynamics of the strands 
at the equilibrium critical temperature T = Tc using Monte Carlo Rouse dynamics. We find that 
the dynamics is anomalous, with a characteristic time scaling as r ~ ^2.26(2)^ exceeding the Rouse 
time ~ L'^'^^. We investigate the probability distribution function, the velocity autocorrelation 
function, the survival probability and boundary behavior of the underlying stochastic process. These 
quantities scale as expected from a fractional Brownian motion with a Hurst exponent H = 0.44(1). 
We discuss similarities and differences with unbiased polymer translocation. 

PACS numbers: 05.40.-a, 82.35.Lr, 87.15. A-, 87.15. H- 



I. INTRODUCTION 

There has been an ongoing interest in recent years in 
the analysis of models of polymer dynamics. The origin 
of this interest is due to two main facts. Firstly, experi- 
ments allow nowadays to control polymers at nanoscales 
and to follow the behavior of single molecules, providing 
thus many insights on their dynamics [IH3] . This has mo- 
tivated more theoretical research in the field. Secondly, 
many of these systems show an anomalous dynamics, a 
paradigmatic example being that of a polymer translo- 
cating through a nanopore |4j. Modeling anomalous dy- 
namics has attracted quite some attention in the Statis- 
tical Physics community due to the ubiquity of this be- 
havior in many physical systems as disordered media [5] , 
conformational fluctuations of proteins [6], diffusion of 
molecules in cells [7] and polymers [B]. 

In the case of polymer translocation, the subdifFusive 
behavior is inferred from the scaling r ^ L" with a > 2, 
which relates the translocation time r to the polymer 
length L. Although there have been a large number 
of publications [H [5HTT] there is no general agreement 
on the value of a obtained from simulations. Also on 
the theoretical side different predictions for a have been 
made [1[I1[T3]. 

Besides translocation, there are other polymer pro- 
cesses which are expected to show anomalous dynam- 
ics. We are interested here in the analysis of zipping 
dynamics, which is the process through which two "com- 
plementary" strands attached by one end bind/unbind 
from each other, as shown in Fig. [l] In our model we do 
not allow bubbles to be formed, hence the dynamics pro- 
ceeds sequentially as in a zipper. We focus here on the 
dynamics at the transition temperature T = Tc- At low 
temperatures (T < Tc) the system is driven towards the 
fully zipped state (Fig. [l|^c)), whereas at high tempera- 
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FIG. 1. (Color online) Sketch of the zipping/unzipping dy- 
namics. The two polymer strands are joined by a common 
end O. Single stranded polymers are shown as dashed lines, 
while double stranded as a thick solid line, (a) Fully unzipped 
state, (b) Partially zipped, (c) Fully zipped. Y denotes the 
position of the end of the double stranded segment. The total 
length of each strand is L. A partially zipped configuration 
is characterized by two single strands of length L — s and a 
double strand of length s. 

tures (T > Tc) the two strands unbind (Fig. [TJa)). The 
dynamics for these two cases were studied in Ref. |14j . 
Interestingly, the zipping time was found to scale as the 
length of the polymer as ~ L" with a « 1.37, whereas 
the unzipping time scaled as Tu ~ L. If the dynamics 
were sufficiently slow (see below), so that zipping and 
unzipping proceed through quasi-equilibrium states, the 
times would scale linearly in the strand length L. This 
corresponds to the motion of a Brownian particle, the 
fork point Y shown in Fig. [ijb), in a linear downhill 
potential. A scaling Tz ^ L" with a ^ 1 implies anoma- 
lous dynamics. We note, in addition, that the anoma- 
lous exponent for zipping (a = 1.37) turns out to be in 
agreement with that found in some simulations of forced 
polymer translocation [13 [H] . 

We recall that in forced (or biased) translocation an 
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external field drives the polymer preferentially towards 
one of the two sides of the separating membrane. This 
can be realized experimentally, for instance for a charged 
polymer as DNA or RNA by imposing an electric poten- 
tial difference on the two sides. In unbiased translocation 
the external field is absent and the polymer translocation 
is driven by thermal fluctuations. 

The aim of this paper is to extend the analysis of the 
zipping dynamics to the critical temperature T = Tc- In 
this case there is no strong bias towards either the zipped 
or unzipped state, as these two have the same free en- 
ergy in the thermodynamic limit. One could therefore 
expect an analogous behavior of that found in unbiased 
polymer translocation. However, the exponent we find 
in our simulations disagrees with those conjectured for 
unbiased translocation suggesting that critical zipping 
is in another universality class. We find however that 
the underlying stochastic process is well-described by a 
fractional Brownian motion [17 for which we determine 
the Hurst exponent, probability distribution function and 
survival probabilities. 



II. MODEL 

The model discussed here was also used in recent stud- 
ies of renaturation dynamics [18] and zipping dynamics 
[H] . Two polymers defined on a face-centered-cubic (fee) 
lattice are joined by a common end. The monomers on 
both strands are labeled by an index i — 0,1,..., L. 
The two strands are self- and mutually avoiding, with 
the exception of monomers with the same index i, which 
are referred to as complementary monomers. Two com- 
plementary monomers can bind by overlapping on the 
same lattice site. A typical simulation run start with 
all monomers < i < L/2 bound, and the monomers 
i > L/2 unbound (see Fig. [l|b)). This initial configu- 
ration is relaxed to equilibrium by means of pivot [19\ 
and local moves which leave the number of bounded 
monomers unchanged. Once an equilibrium configura- 
tion is obtained the actual simulation is started. The 
pivot algorithm is no longer used and the Monte Carlo 
updates are strictly local. They consist of corner-flip 
or end-flip moves that do not violate self- and mutual 
avoidance. A feature of the fee lattice geometry 
is that the local moves are such that two differ- 
ent polymer chains cannot intersect each other. 
In the simple cubic geometry the non-crossing 
condition can be realized by the bond-fluctuation 
model [20j. The local dynamics on the fee lattice 
reproduces the Rouse model behavior, as shown 
in [21\. 

A Boltzmann weight a; > 1 is associated to the binding 
of two monomers. In the Monte Carlo dynamics binding 
occurs with probability 1 while unbinding with proba- 
bility 1/uj, so that detailed balance is satisfied. A Monte 
Carlo step consists in selecting a random monomer on one 
of the two strands. If the selected monomer is unbound a 



local fiip move is attempted. If the selected monomer is a 
bound monomer there are two possibilities. Either a local 
fiip of the chosen monomer is attempted, and if accepted, 
this move results in the bond breakage; or a fiip move of 
both bound monomers is generated, which does not break 
the bond between them. In the model discussed here we 
do not allow any bubble formation neither for zipping 
nor unzipping, by imposing the constraint that monomer 
i can bind to its complement only if monomer z — 1 is 
already bound. Analogously monomer i can unbind only 
if monomers i + l are already unbound. This is the model 
Y which was referred to in Ref. [H]. As unit of time we 
take 2L + 1 Monte Carlo steps, so that in one time step 
one attempted update per monomer is performed. 

III. EQUILIBRIUM FREE ENERGIES 

We start by discussing some properties of the equi- 
librium free energy of zipping polymers. We consider in 
particular the dependence of this free energy on the coor- 
dinate s, describing the position of a fictitious Brownian 
particle. We analyze the motion of this particle in the 
given free energy landscape. 

The number of configurations for a linear polymer of 
length L, in the limit of i — > oo, takes the asymptotic 
form Z ^ where /i is the connectivity constant 

and CTi a universal exponent associated to the end ver- 
tex of the polymer [2? (we use here a different notation 
from the customary exponent 7 = 1-1- 2ai). In gen- 
eral the partition function of polymer networks of more 
complex topology [33] is also characterized by sublead- 
ing universal exponents. For instance a star polymer with 
three arms of length L has a partition function scaling 
as Z3 ~ /i'^^L'^^+^°'i , with 0-3 the exponent associated to 
a vertex with three outgoing arms. Here the factor Scti 
accounts for the three end vertices. 

We can now consider the case of partially zipped poly- 
mer strands of Fig. [l|b), which is related to that of a 
three arms star polymer. We allow now the three arms 
to have lengths L — s and s. In addition we have to 
account for the Boltzmann weight uj associated to the 
binding between the s monomers on the double stranded 
segment. The total partition function is then given by 

Z{s) = fi^^-'uj'L^'''+'''f{s/L) (1) 

where f{x) a scaling function. We analyze now the lim- 
iting scaling behavior of the partition function Q. For 
s — > one should recover the partition function of a sin- 
gle polymer of length 2L which yields f{x) ~ x"^'^'^^ . 
The analysis of the limit s ^ L similarly imposes f{x) ~ 
(1 — xY^'^'^^ . One can incorporate the two limits in the 
following expression: 

(2) 

with g{x) an analytic function of its variable. 
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We distinguish now two cases: uj ^ [i. and uj = fj,, 
which correspond to the off-critical and to the critical 
case, respectively. If lo ^ fi the leading dependence on s 
of the free energy is 



f{s) = - log Z{s) slog(a;/^) 



(3) 



where the constant terms, and subleading terms in s have 
been omitted. For lo > fj, the zipped state (low temper- 
ature phase) is favored, whereas lu < ^j, favors unzipping 
(high temperatures). Viewing s as a zipping coordinate 
one can consider a Fokker-Planck (FP) description of the 
process, which is the motion of a Brownian particle on 
a potential /(s). For uj ^ fj, the potential is linear in 
s (Eq. ([3])), equivalent to a biased Brownian motion to- 
wards s — (uj < fi) or s — L {io > In both cases 
the characteristic time of the particle for starting in one 
boundary to reach the opposite one scales as ~ L. 
At the critical point (a; = /i) the free energy becomes 



f{s) = -((Ti +0-3) log 



1 - 



L 



- log gis/L) (4) 



(omitting terms which do not depend on s). We compare 
now this expression with the free energy of translocating 
polymers. Consider a polymer on a pore of an infinitely 
wide separating plane. The pore divides the polymer 
into two non interacting parts of lengths s and L ~ s 
respectively. The two parts are characterized by parti- 
tion functions Zi ^ fj,''s'^''~^ and Z2 ^ ii'"^^{L — s)^=~^, 
where 7^ is an entropic exponent associated to a polymer 
attached to a planar surface. One obtains the equilibrium 
free energy of the polymer on the pore ^ : 



/t(s) = -(7.-1) log [^( 



1 - 



(5) 



This free energy has a strong analogies with that of the 
zipping polymer at Tc (Eq. |4]) . One difference is the value 
of the exponent. For a three dimensional self-avoiding 
walk attached to a planar surface 7s = 0.70 [55], whereas 
from rcnormalization group results |24j one finds ui + 
CT3 « -0.08. 

In an early study of translocation, Chuang et al. [3] 
analyzed the FP equation corresponding to a Brownian 
particle moving in the potential ([5]). Within this frame- 
work they showed that, by appropriate rescaling of poly- 
mer length L, time and space coordinates one finds a 
translocation time scaling as rf^ ~ L^, independently 
on the value of 7^. A similar rescaling is also applicable 
to the free energy in Eq. (|4|, even in the presence of an 
analytic scaling function g{x). The analysis thus shows 
that the zipping/unzipping time 
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L . However it 

was also argued [3] that this result is not self-consistent 
as the time required for a polymer to equilibrate (the 
Rouse time, ~ L}^"^^ and v « 0.59 in three dimen- 
sions thus Tji ^ l?-^^) turns out to be longer than the 
predicted ^ from the FP equation. In conclusion 
a polymer cannot translocate through quasi-cquilibrium 
states, because it would do so at a rate at which it can- 
not equilibrate A similar reasoning can be applied 
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FIG. 2. (Color online) Log-log plot of the zipping (r^) and 
unzipping (r„) times as functions of the strand lengths L. 
The two quantities scale with the same anomalous exponent 
Tz ~ Tu ~ ^2.26(2) Qij-^igg ^jj^j squares are data from simu- 
lations with initial condition sq = L/2, i.e. half-zipped con- 
figurations. Diamonds correspond to the zipping time using 
as initial condition so — 3L/4. Triangles correspond to the 
unzipping time with initial condition so = L/A. 



to the potential of zipping strands (Eq. Q). It can be 
shown that the modulating scaling function g(x) does 
not alter the general result and thus the FP equation 
predicts a zipping time scaling as r ^ . As for translo- 
cation, however, this result is not self-consistent since the 
equilibration time of the polymer would then exceed the 
zipping/unzipping time. One deduces also for zipping at 
Tc a lower bound t^/u ^ L^'^^'^ of the zipping/unzipping 
time. 



IV. SIMULATION RESULTS 
A. Scaling of zipping/unzipping times 

The simulations at T — Tc are started from a configu- 
ration with s = So, i.e. with a star polymer with a dou- 
ble stranded segment of length sq and two single strands 
of length s — So- The critical point is obtained by set- 
ting the Boltzmann weight of binding of two monomers 
to w = /Lt. The connectivity constant for self-avoiding 
polymers in a fee lattice is known with a high degree of 
accuracy /i — 10.0362(6) PP, which provides T^. with a 
very good precision. The initial configuration is thermal- 
ized while keeping s fixed. This constraint is released 
at time t — and the simulation is stopped once one of 
the two boundaries s = or s = L is reached. These 
boundaries correspond to fully unzipped strands and to 
fully zipped strands, respectively (see Fig. [l]^a) and (c)). 
The simulations are repeated typically for 10^ different 
realizations. 

Figure [2] shows a plot of the averages of the zipping 
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and unzipping t„ times as a function of the strand length. 
Error bars are smaller than symbol sizes. The four dif- 
ferent data sets correspond to different initial conditions 
So = i/2 (half-unzipped), sq — L/4 and sq — 3L/4. 
Apart from small deviations, for very short polymers, 
the data show a power-law behavior which is described 
by the same exponent ^ t„ ~ ^2.26(2)^ Note that 
for So = L/2 the unzipping is systematically slightly 
slower than zipping (tu > t^). The differences are due to 
the asymmetry of the problem with respect to the inter- 
change s — > L — s. This asymmetry has seemingly not a 
strong effect on the dynamics; it only influences and 
T„ by a multiplicative factor, but it does not change the 
exponent. The exponent a = 2.26(2) indicates a subd- 
iffusive behavior. In addition the expected lower bound 
Tji ^ L^+^'', discussed above, is verified. But clearly 
a > 1 + 2v ^ 2.18. Also notice that a is considerably 
lower than the value 2 -\- v = 2.58 recently conjectured 
for unbiased translocation [T^ . 



B. Probability distribution function (pdf) 

We computed next P(s, t), the probability distribution 
function (pdf) of the value of s at time t. Initially, at time 
t = 0, one has -P(s, 0) = 5{s — so), whereafter P{s,t) 
spreads in time. Here we consider the case sq = L/2. 
The shape of the distribution at later times gives some 
insights on the nature of the underlying stochastic pro- 
cess. It has been recently suggested [26^^ that unbiased 
translocation could be described by a class of processes 
known as fractional Brownian motion (fBm). The fBm 
|17j is a Gaussian, self-affine process with stationary in- 
crements. It is described by a probability distribution 
with a Gaussian shape, but with variance growing as a 
power law in time, i.e. 



P{s,t) 



1 



y/2-KDt^" 



exp 



(6) 



where so = L/2 in our setup (the starting point of the 
simulations) and is a constant. Here < iJ < 1 is 
known as the Hurst exponent; H = 1/2 corresponds to 
Brownian motion. H < 1/2 and H > 1/2 are the subdif- 
fusive and the superdiffusive cases. The distribution ([6| 
cannot hold at all times, due to the presence of bound- 
aries at s = and s = L. 

Figure [3] (main graph) shows plots of P(s, t) /P{L/2, t) 
vs. s on a log-linear scale at different times for strands of 
length L = 192, averaged over 10^ histories. In these sim- 
ulations, P(s, t) is calculated from the surviving samples, 
i.e. those which never reach the boundaries s = and 
s = L at time t. Similar results were obtained for other 
L-values. One observes the spreading of the distribution 
from its initial delta shape. For large times, P{s,t) con- 
verges to a stationary distribution whose properties we 
will discuss below. 

As a check whether (and when) these data can be de- 
scribed by a Gaussian distribution we have calculated the 




FIG. 3. (Color online) Plot of P{s,t)/P{L/2,t) for strands 
of length L — 192 for different simulation times t = 2^, 2^^, 
2", 2^^, 2l^ 2", 2^^ and 2^*. The inset shows P{s = L/2,t) 
versus time in a log-log plot. 
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FIG. 4. (Color online) Plot of 7? = {(s - (s))*)/((s - {s)ff 
as a function of time for various L- values. The dashed line 
indicates the Gaussian value R = 3. The sampling is over 
5 ■ 10^ configurations for L = 24, 48, 2 • 10^ configurations for 
i = 96 and 10^ configurations for L = 192, 384, 768. 



ratio R= ((s— (s))'*)/((s— (s))^)^ which for a Gaussian 
equals 3. Results for i? as a function of time and for var- 
ious values of L are shown in Figure |4j These data show 
that after an initial regime the gaussian value is very 
slowly approached. After a time of the order of the zip- 
ping/unzipping time, the boundaries are being felt, and 
a deviation from the Gaussian value develops, as should 
be expected. The extrapolation for different values of L 
yields R = 3.04(5), consistent with a gaussian distribu- 
tion. 

We next turn to the determination of the exponent 
H. This can be done in two ways. Firstly, we looked 
at P{s = L/2,t) which should decay as . As can be 
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FIG. 5. (Color online) Plot of as a function of time t for 
different values of L. The dashed line is a power-law fit to the 
data, yielding H = 0.44. 




FIG. 6. (Color online) Log-log plot of minus the velocity 
autocorrelation as a function of time. The dashed line is a 
straight line (not a fit) with slope —1.12. 



seen in the inset of Fig. [3] our results are well fitted with 
a power law. From an analysis of the data for various 
L- values we arrive at the estimate H = 0.44(1). 

Secondly, we determined the scaling of the variance of 
the position distribution, o'^(i) = ((s — (s))'^), as a func- 
tion of time. For the Gaussian distribution (|6| one has 
cr^(i) ~ t^^ . For times above Tz/u deviations from this 
behavior appear and asymptotically in time should 
saturate at a value ^ L^. We therefore expect a scaling 
of the form 



a' it) - t^^Fit/r) 



(7) 



where F{x) is a scaling function and r is or r^. 

This is precisely the behavior that we observe in our 
simulations (Fig.jsj). From the initial power law increase 
we arrive at the independent estimate oi H ^ 0.44(1). 
Since at i w r, cr^ ~ we obtain ~ r„ ~ where 
a = 1/H. For H = 0.44(1), a = 2.27(2), fully consistent 
with our results on the behavior of the zipping/unzipping 
times shown in Fig. [2] In conclusion, we find clear numer- 
ical evidence that the dynamics of the zipping coordinate 
s is well described by the distribution ^ with an expo- 
nent H = 0.44(1) in the time regime where t < t^/^. 



C. Velocity autocorrelation 

Since fBm is a Gaussian process, it is fully character- 
isted by its average and covariance matrix. The latter, 
which can also be interpreted as a position correlation 
function, has the form 17 



{siti)sih)) = Diti" + 4" - \h - t2r) (8) 

where D is a constant. The autocorrelation of the ve- 
locity V — ds/dt can be obtained by differentiation with 
respect to ii and t2 [28] 



{Vih)vit2)) 



dtidt2 



d'^\ti-t2 
dtidt2 



\2H 



2H{2H - l)D\ti -i2p^"^ 



+2HD\h-t2\'"-'S{h-t2) 

Because of the stationarity of the process one obtains for 
t > 



{v{t)v{Q)) = 2H{2H- l)Dt 



2H-2 



(9) 



For a fBm with H < 1/2 (subdiffusive), as is the case for 
zipping, this implies a negative velocity autocorrelation 
function. Using H = 0.44 one finds 2H -2 = -1.12. In 
Fig. [6] we show our results for this quantity in a log- log 
plot. In order to compute (u(t)u(O)) reliably more than 
10^ configurations were sampled, which is 100 more than 
typically done for other quantities. Still, the velocity au- 
tocorrelation could be estimated with sufficient precision 
only for short times [t < 10^). It is therefore difficult to 
obtain a good estimate of H from these data. The fig- 
ure shows however that {v(t)v(0)) is negative and decays 
with a power law as expected. The decay for i > 10 is 
not inconsistent with the predicted value —1.12. 



D. Survival probability 

Next, we look at the probability distribution of the 
(un)zipping time Q^r^/y^, L). These results are plotted 
in Figure [t] for L = 24, 48, 96 and 192 (only unzipping 
data are shown, the distributions for the zipping times 
are very similar). As seen from the figure, the proba- 
bility distribution decays exponentially for long times. 
The data for different lengths collapse into a single scal- 
ing function (for L sufficiently large), where t/L?'"^^ is 
used as a scaling variable. The normalization condition 
of the probability distribution function implies the fol- 
lowing scaling form: 



Q{T,L)^L-'-'^f{TlL 



-2.26\ 



(10) 
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where f{x) is a scaling function which decays exponen- 
tially for large x. An exponential decay of the sur- 
vival probability S{t), which is related to the pdf of the 
(un) zipping time by Q(t) — —dS/dt is also characteristic 
for fBm with two absorbing boundaries . These re- 
sults also rule out a description of P{s,t) in terms of a 
fractional Fokker-Planck equation with absorbing bound- 
aries for which Q(r) decays as a power law [29. .30] , 



E. Boundary exponents 

There has been quite some interest in the behavior 
of fractional Brownian motion in the vicinity of an ab- 
sorbing boundary. For ordinary diffusion the probability 
distribution P{s,t) on the interval [0,L] becomes pro- 
portional to sin(7rs/L) asymptotically. This means that 
P(s, t) vanishes linearly with the distance y from an ab- 
sorbing boundary. For a fBm with Hurst exponent H 
it has been recently argued [26] that the pdf vanishes 
as ^ y'^, with (j) = {1 — H)/H. For ordinary diffusion 
{H — 1/2) this gives the correct result 0=1. 

In view of this recent interest, we have investigated 
the behavior of the pdf near absorbing boundaries. In 
order to determine the boundary exponent we have 
replotted in Figure [8] the results of Figure [3] in a log- 
log plot. According to Ref. [26] a fBm with a Hurst 
exponent H — 0.44 would imply a boundary exponent 
(j) = 1.27. The corresponding power-law behavior is 
shown as a dashed line in Figure [H] This exponent is in 
good agreement with the numerical data at sufRciently 
long simulation times for which the distribution becomes 
stationary. 
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FIG. 7. (Color online) Linear-log plot of the scaling function 
of the probability distribution of the unzipping time for the 
sizes L = 24, 48, 96 and 192. The calculations are done 
over 10^ configurations for L =24, 48 and 96 and over 2 ■ 10^ 
configurations for L = 192. 
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FIG. 8. (Color online) Plot of InP(s,t) vs. Ins emphasizing 
the vanishing of the position pdf in the vicinity of the absorb- 
ing boundaries (same data as in Fig. [3) L = 192 and times 
are t = 2^, 2^°, . . . , 2^"^ , 2^*). For the two longest simulation 
times the distribution is almost stationary and the probabihty 
distribution vanishes in agreement with a power-law behavior 
~ s^'^^ (dashed line). This is the expected behavior for a 
fractional Brownian motion with H = 0.44 in the vicinity of 
a boundary, as derived in Ref. |26j . 



V. CONCLUSIONS 

In this paper we have analyzed the dynamics of zipping 
polymers at their critical point by means of simulations 
of polymers undergoing Rouse dynamics. The natural co- 
ordinate describing the process is s, denoting the length 
of the double stranded (zipped) part of the strands. At 
T — Tc the equilibrium free energy has a form very sim- 
ilar to that of a polymer translocating through a mem- 
brane, i.e. a logarithmic dependence on s, except for a 
weakly modulating function, which is shown to have a 
weak effect on the dynamics. 

Using the same arguments as for translocation one can 
derive a lower bound r^/^ ^ L" with a > 1 -I- 2j/ w 2.18. 
The numerical results are indeed in agreement with this 
behavior, and provide an estimate a sa 2.26(2) which is 
higher than the lower bound value. Different exponents 
have been reported in the literature for polymer translo- 
cation (ranging from 2.23 [S], 2.51[21] to 2.59 [H]), but 
it should be noticed that critical zipping does not seem 
to share a common universal dynamical scaling behavior 
with unforced polymer translocation. This could be a bit 
surprising at first, as a recent study of zipping dynamics 
at T < Tc showed instead a good agreement with that 
found in forced polymer translocation [T3] . On the other 
hand it also reported that that the unzipping dynamics 
at T > Tc is not anomalous (t„ ~ L). At T = Tc the 
dynamics is a complex combination of zipping and un- 
zipping resulting in an exponent that deviates from that 
for unbiased translocation. 

It has been recently suggested that translocation could 
be described by fractional Brownian motion ^25] [571 [33] . 
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We have therefore considered the possibihty that a fBm 
approach could also provide a consistent description of 
the critical zipping/unzipping dynamics. Our results 
show that this is indeed the case. On time scales smaller 
than the zipping/unzipping time, the distribution of the 
zipping coordinate s is well described by a Gaussian with 
a variance that grows as with H = 0.44(1). In the 
presence of absorbing boundary conditions we find that 
P{s, t) conditionned on not being absorbed yet, converges 
to a distribution that behaves as a power law near the 
boundaries. The associated exponent agrees with that 
recently determined for fBm [26] . 

Fractional Brownian motion has been observed in other 



polymer processes, most notably in the fluctuations of the 
distance between an electron transfer donor and accep- 
tor pair in single molecule experiments on proteins [5]. 
Recently, it has been argued that this complicated dy- 
namics can arise from a superposition of Markovian fluc- 
tuations of normal modes as for example those present 
in the Rouse model fB?. It would be interesting to get 
further insight in the underlying dynamics of the zipping 
coordinate s in a similar spirit. 
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